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We show that the evolution of magnetic fields in a primordial plasma, filled with Standard Model 
particles, at temperatures T > 1 OMeV is strongly affected by the quantum chiral anomaly - an effect 
that has been neglected previously. Although reactions equilibrating left and right-chiral electrons 
are in deep thermal equilibrium for T < 80 TeV, an asymmetry between these particle develops in 
the presence of strong magnetic fields. This results in magnetic helicity transfer from shorter to 
longer scales. This also leads to an effective generation of lepton asymmetry that may survive in 
the plasma down to temperatures T ~ f MeV, which may strongly affect many processes in the 
early Universe. Although we report our results for the Standard Model, they are likely to play an 
important role also in its extensions. 
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Magnetic fields are expected to play an important role 
in the early Universe Recent observational indications 
of the presence of magnetic fields in the inter-galactic 
medium [IHII suggest that cosmological magnetic fields 
(CMF) may survive even till the present epoch. Thus 
they could have played the role of seeds for the forma- 
tion of galactic magnetic fields. A number of mechanisms 
for the creation of CMF at very high temperatures have 
been proposed (see e.g. 041] and refs. therein). 

In this paper we concentrate, however, on a different 
problem: we assume that strong CMF were already gen- 
erated at a temperature > 100 GeV and we study the 
subsequent evolution of such fields. Usually, this evolu- 
tion is described by the system of Maxwell plus Navier- 
Stokes equations (for a detailed review see @, 0)- Here 
we will argue that, for temperatures T > 10 MeV, this 
system of MHD equations should be extended to include 
a new effective degree of freedom, even if all particles and 
reactions arc described by just the Standard Model of 
particle physics. This significantly affects the evolution 
of CMF and the state of the primordial plasma. 

At such temperature rates, of all perturbative pro- 
cesses related to the electron's finite mass are suppressed 
as (m e /T) 2 . Ignoring these corrections for a moment, the 
number of left and right-chiral electrons 1 is conserved in- 
dependently. 2 That is, apart from the vector current 
= i/i^tp describing conservation of electric charge 
(?il + n R ), the average number density of the left- (right) 
chiral electrons jil,r = I d 3 xip <! (l ± 7s)V> does not 



1 More precisely, is the difference between the number of left 
particles and left anti-particles (same for nji). 

2 The number of left-chiral electrons is not conserved when weak 
processes are fast (the conserved quantities are n/, + n„ e and 
nn). The coefficient in Eq. below takes this into account. 



change with time. This is true on time scales smaller than 
the chirality- flipping scale TJ 1 . Although the chirality- 
flipping rate is suppressed as compared to the rate of 
chirality-preserving weak and electromagnetic processes, 
it is faster than the Hubble expansion rate, H(T), for 
temperatures below 80 TeV Q and chirality nipping pro- 
cesses are in thermodynamic equilibrium. Yet on time 



scales T, 



< t < r. one should introduce inde- 



EM,weak ^ " ^ x / 

pendent chemical potentials, (il and fx R , for two approxi- 
mately conserved number densities, with til,r = ^'T 2 . 
In the presence of external classical fields the conserva- 
tion of the axial current is spoiled, however, by the chiral 
anomaly Q - a quantum effect leading to a change of 
n L - n R : 



d(n L - n R ) 
dt 
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where a = 4- is the fine-structure constant and rt is the 
magnetic helicity defined as 

W(t) = i J y d3xA - B . ( 2 ) 

(where B is the magnetic field and A the vector potential, 
with B = V x A). The quantity ^ is gauge invariant, 
provided that B is parallel to the boundary of V (sec 
e.g. 0)- The time evolution of H(t) is given by Q 



dU 
dt 



2 
V 



d 6 xE-B . 



(3) 



In terms of the difference of left and right chemical po- 
tentials, Afi = fi L — fi R , Eq. JTJ reads 



d(Afi) c A a dU(t) 



(4) 



dt T 2 dt ' 

where ca is a numerical coefficient of order one that 
describes the dependence of nr, on globally conserved 
charges in the primordial plasma. 
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If Ap 0, the chiral anomaly leads to an additional 
contribution to the current in Maxwell's equations (Tol — 
0: 



V x B = aE + -Ap(t)B 

TT 



or, combining it with the Bianchi identity Vx£ = — B: 



at a 



(5) 



As weak reactions are fast enough at these temperatures 
to establish local thermodynamic equilibrium (LTE), 
in the background of long-wavelength electromagnetic 
fields, space-dependent chemical potentials pl,r(x) may 
be defined. Eqs. (fTJ), (Q| can then be written in a local 
form, and Eq. <j5j> acquires additional terms, proportional 
to the gradients of Ap(x) [II [II [II]. We assume fields to 
be slowly varying and neglect these effects as well as those 
depending on the velocity field. We will show that, even 
in this limit, the evolution of magnetic fields significantly 
changes as compared to the usual Maxwell equations. A 
more realistic analysis should include all the derivative 
terms, as well as the Navier-Stokes equation describing, 
in particular, turbulent effects known to be important for 
the evolution of CMF. We leave a more complete micro- 
scopic derivation and an analysis of the full system to 
future work and use the simple model described above to 
illustrate the previously neglected effects. 

Eqs. ((3HS1) remain valid in an expanding Universe if 
written in conformal coordinates (see e.g. 0, [U, fig)- 
Henceforth we use conformal quantities and define con- 
formal time as r\ 



4p, where M* 



is the effective number of relativistic degrees of freedom. 

Eqs. (J4HSJ) are translation and rotation invariant. We 
introduce the magnetic helicity density, Hk , and the mag- 
netic energy density, pk, in Fourier-space, with pb(v) = 
J dk pk(ri) and "H(r/) = J dkHkiv)- 3 The quantities Hk 
and pk obey the inequality \Hk\ < jPk, which is satu- 
rated for field configurations known as maximally helical 
fields. In our subsequent analysis, we focus on this case 
and choose for definiteness Hk > and Ap > 0. Mul- 
tiplying the Fourier version of Eq. ([5]) by the complex- 
conjugate mode B^, we obtain, after some simple manip- 
ulations (see Appendix [El for details, cf. [H|), 



dH k 

Ot] 

d{Ap) 
dr/ 



2k 2 a kAp 
H k + -Hk , 

0~ c 7T 0~ c 

-(c A a) Jdk^-^-TfAp 



(6) 



where we have restored the chirality flipping rate T f in 
Eq. and used the conductivity er c = cr(?y)/T « 70 (2lj. 

The system ([6HZ]) has been previously studied in two 
regimes. It was demonstrated in [HI, dHHH that, in the 
presence of a large initial chemical potential difference 
Ap(rj) > 0, the quantity 



74(77) =?^exp{^(|- J\(fj)dfj - k(r, 

(8) 

grows exponentially fast for sufficiently long wavelengths. 
Conversely, in [lU the initial background of helical (hy- 
per)magnetic fields was used to generate a non-zero 
chemical potential for T > 100 GeV. 

In this work, however, we consider helical CMF with 
some initial spectrum Tifj,, already present at T ~ 
100 GeV in the hot plasma, filled with particles in ther- 
mal equilibrium (cf. [l8l - l20l . |23| ) . It was believed that as 
T f > H(T) for T < 80 TeV, no chiral asymmetry will 
survive. 

Chirality evolution. Below we show that both Ap 
and the magnetic helicity do survive below 100 GeV on 
time-scales much longer than diffusion or chirality flip- 
ping times (till 10 + 100 MeV). (For T f -» the sys- 
tem (|BH7|) can even reach a stationary state with non-zero 
B and Ap). 

To see this, it is convenient to separate on the right 
hand side of Eq. ([7]) a source term 5b (77) (independent 
of Ap) (cf. E3): 



Mpi and 3* where 



d(Ap) 
dr] 

r B (r?) = 

Sb(v) = 



-(T B (r,) + T f )Ap + S B (r,) 



(9) 
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dk k%k = 



dkk 2 n k 



2ca« 2 
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We begin our analysis of Eqs. and (| MTU)l with the 
case where = and the field is initially "monochro- 
matic", i.e., 



H k (r 1 ) = n(v)S(k-k ) 



(11) 



The form (fTTj) is preserved during the evolution as Eq. (|6j 
is homogeneous. 4 Putting in Eq. (|9]) d(Ap)/dr] = and 
Tf = we find the so-called tracking solution: 



Apt 



5b(t7) _ 27rfc 



(12) 



This is an exact static solution of the system (O, ©: 
A/itr arid H(ri) remain constant, i.e., dissipation due to 



In terms of left and right circular polarized modes 73 fc , = 
2^(|B,t| 2 ~ \ B k\ 2 ) Pk = (^t(|B+| 2 + |73-| 2 ). Integrals 
over k run over the radial direction only (see e.g. llgU20lO . 



4 This is an artifact of our homogeneous approximation with 
Afi independent on spatial coordinates, see below. 
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magnetic diffusion is exactly compensated by growth due 
to a non- vanishing chemical potential difference A/i tr ; 
(cf. ©). 

Until now we have completely neglected the massive- 
ness of the electrons. It is straightforward to compute 
that the rate r /(r?) due to electromagnetic processes is: 5 

17(77) w a 2 (jjfc) V 2 - E q s - ©, ® can be rewritten 
to describe deviations from the equilibrium static solu- 
tion (IT21): 



dAp 



di] 

dT B _ T_b_ 
di] n a 

_ 2A: 2 



r s (A/i-A^ tr )-r / A^, 

A/j 



A/J tr 



1 



(13) 
(14) 



where rj a = is the magnetic diffusion time. From 
Eq. (fT3")) we see that Tb and T/, which enter symmetri- 
cally in Eq. (|9]), play very different roles. The rate Tf, 
that depends only on temperature, constantly drives Ap 
to zero. The Tb term pushes the system towards the 
equilibrium value (| 12|) (that depends only on ko). It de- 
pends on pb and has its own dynamics, Eq. (|14p. 

If the magnetic field is large (such that Tb 3> 17), 
any initial value of A/z will be quickly "forgotten" and 
Ap will be driven towards Apt r - At that moment a new 
tracking solution will take over, with Ap— Ap.tr ~ jAfM, T , 
where 



r s (r?) 



(15) 



This new solution is valid, provided two conditions hold: 
(%) 7 <C r B i and (m) 7 <C ^bVct- When this holds the 
evolution of Tb is given by (fl~4")) . 



d?7 



1 

la 



(16) 



Eq. (|16[) shows that remains practically constant 
when 77 <C 7 loll{ 7 i)i which is significantly longer than 
r\ a . as 7 -c 1. To estimate the time at which the function 
7(77) ~ 1 we note that it evolves with time because of an 
increasing chirality flipping rate T f (rf) oc ?7 2 and because 
the total magnetic energy dissipates (fT!))) . Neglecting this 
latter change, we estimate 7 to be given by: 



7TO- c 

2c A 



3MJ «ig,r B 



1CT 5 ZlOOMeV 



re 



(17) 

where we used tb = pb/(t$q9*T ) is the fraction of mag- 
netic energy density to the total energy density. From 



Eq. (|16[) we see in addition that Tb remains approxi- 
mately constant as long as ^- J j(n)drj < 1. Using (fT?) 
we find that (this is illustrated in Fig. lAll in Appendix[A| . 



377 CT 



< 1 . 



(18) 



Inverse cascade. So far we have considered a toy model 
example of " monochromatic" helical field (fTT|) . Although 
the Eq. ((5]) is linear, the modes Hk are not independent 
for different k (due to the integral in the Eq. (J7J). For 
a continuos spectrum, this interaction results in another 
very important effect: the initial spectrum reddens with 
time, the total helicity being conserved (similarly to the 
"inverse cascade" phenomenon Q, Sec. 7.2.3. 

Indeed, let us consider first the case of two modes 
(fci, "Hi(?7)) and {k 2 ,H 2 (rf)) with k\ > k 2 , to understand 
the situation qualitatively While Tb Tf, the evolu- 
tion for Ap has the form: 



d{Ap) 
drj 



kiUi+k2U 2 )Ap l +^^(kl'H 



-k\Ui 
. (19) 

One can again try to construct a tracking solution of 
Eq. (fT9|) by putting its l.h.s. to zero. It is clear, however, 
that, unlike in the case (fl2"j) . such a tracking solution 
cannot be time independent. Indeed, according to Eq. ([5]) 
Hk = only if Ap = while our solution A^ t r = 

— i^TT-TT^Tr depends on both modes. In the case where 
a shorter mode (fci) contains most of the energy density, 
Ap will grow very fast and reach: Ap^ rj ?lEJSx{\ _ e ) 
Initially e = f 2 ^ 2 < 1 as H 2 is subdominant and Ap is 
close to its "static" value for k\ . Therefore the mode Hi 
remains almost constant, for 77 < rj a {ki) / e{rj). For the 
mode H2, however, k 2 < a ^ tr and from the solution ([§]) 
(valid for any Ap(j])) we see that H 2 will start growing. 
As its growth enters the exponential phase, e increases 
and k\ becomes greater than Ap(r\), causing Hi to decay 
exponentially Ap(n) will therefore quickly evolve to the 
value ^fc 2 . From Eq. © we find that for e < 1: 



(20) 



the 



and see that Hi = —H 2 as long as Tb 3> Tf, i.e. 
total helicity of the system is conserved. 6 

The evolution of continuous spectra is qualitatively 
very similar. Assume that the initial helicity spectrum 
H\ has its maximum at a scale k\ and then decays as 
Hi oc (t^)™ s ~ 2 , with n s > 3. The scale ki determines 
the value of Ap at the beginning, while the longer modes 



6 H is related to A/i by Eq. J1J. Even if A/i changes smoothly, the 
5 The contribution of weak processes to Tf is small at T < very small numerical coefficient in J4f suppresses the change of 

100 GeV, see Appendix [D] H. 
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FIG. 1: Evolution in the absence of chirality flip for continuous spectrum oc fc 3 . Left: the evolution of A/i(rj) (black solid 
line). The red dashed line: the approximation of the spectrum by 10 modes (k n = 10 -10 x 2"~ 1//2 , n = 1, 10). The horizontal 
lines are tracking solutions for individual modes k n . Right: transfer of helicity from the shorter to the longer modes (red to 
blue). 



grow. At the moment when back-reaction of these grow- 
ing modes on A/z becomes non-negligible, the chemi- 
cal potential difference gets smaller and the modes with 



ki > k > 



2ir 



start decaying. 



For discrete spectrum A/i changes by "steps", defined 
by the modes k n (see Fig. Q] for Tf = 0, red dashed 
curve). Every such "step" corresponds, on the other 
panel of Fig. [1] to a fast decay of one helicity mode and 
exponential growth of an adjacent one (from red to pur- 
ple), while the total helicity remains constant (black dot- 
dashed line). The conservation of helicity implies that 
the total magnetic energy gets dissipated as = ^Hk 
for helical fields. If we sample the same spectrum with 
a larger number of modes, the evolution of A/i becomes 
monotone. If the initial spectrum is sharp (n s > 3), A/i 
will decay more slowly, and the short modes will survive 
for a longer time. The resulting spectrum will, however, 
be roughly the same - the helicity concentrates around 
the longest mode ki that had enough time to start grow- 
ing, A/i = A/i tr (/c2), the magnetic energy is smaller by 
the factor I 2 -. We believe that these results correctly 
describe the interaction between different helicity modes 
even in the inhomogeneous case, provided that deviations 
from LTE are not very dramatic and Afi(x) is smooth. 

Finally, the exact numerical solution of the full sys- 
tem with continuous spectrum and finite 7(7/) is shown in 
Fig. [2] where the red dashed line shows A/j,(j]) in the case 
Tf = and the thick blue and green lines show A/i(r/) 
for Tf 7^ and different r B - This full evolution follows 
that of r f = and then breaks down exponentially fast 
when 7(77) ~ 1 and ([18]) holds. 



Conclusion. This work demonstrates that the tradi- 
tional MHD equations should be modified, when applied 
to a plasma of relativistic particles with T 3> m. The 
proper account of the chiral anomaly changes the evo- 
lution of magnetic fields in two ways: (i) the magnetic 
fields survive several orders of magnitude longer than the 
time defined by magnetic diffusion (Eqs. (jTij)) - ([IT]) ) and 




100 50 



FIG. 2: Evolution of the chemical potential for tb = 5 x 10 5 
(solid blue line). Vertical line at T ~ 150 MeV marks 7 = 
1, A/i(150 MeV) ~ 2.9 x 10~ 6 . The red dashed line shows 
evolution of the chemical potential for Tf =0. The green 
short-dashed line shows A/i(r/) for r B = 5 x 10~ 4 . The purple 
dashed-dotted line shows the decay of the chemical potential 
in the absence of magnetic field.) 



(ii) an "inverse cascade" develops, transferring energy 
from shorter to longer wavelength modes. The effect 
depends on the energy of magnetic fields parametrized 
by its ratio to the total energy density, r B . In the lit- 
erature discussing the evolution of magnetic fields (see 
e -g- [3 US EH) r B < 1 is often considered. It was 
demonstrated e.g. in Refs. (H,^ that r B ~ few x 10~ 3 
may be generated at cosmological first order phase transi- 
tions. The mechanism of [26[ predicts maximally helical 
magnetic fields with B ~ 100 GeV 2 (i.e. r B ~ 10~ 2 , 
c.f. [13) at small scales. See [|, MM, HI HI and refs. 
therein. 

We refer to the value of r B at T - 100 GeV. Due to 
the inverse cascade and helicity conservation this energy 
decreases by about an order of magnitude by the time 
when our effect stops (T ~ 10 - 100 MeV). The sub- 
sequent evolution of the magnetic fields is described by 
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the conventional MHD 0, [IH, a significant part of 
rg further dissipates, so that only the large scale tail of 
the spectrum may survive due to turbulent effects. To 
predict the final fate of these CMF for every initial spec- 
trum and compare them with cosmological bounds (see 
e.g. our results should be combined with the MHD 

analysis. Nevertheless, the above-described mechanism, 
based entirely on the Standard Model, clearly improves 
the chances of survival of CMF generated at subhorizon 
scales HH- Indeed, even for r B ~ 10~ 5 the fields survive 
down to T ~ 100 MeV, while for rg ~ 0.1 the inverse 
cascade is operational down to T ~ 10 MeV. Moreover, 



regardless of the survival of the CMF, this effect is im- 
portant as the left-right asymmetry in the electron sector 
survives down to T ~ 0(100 ) MeV and thus potentially 
affects important processes in the early Universe: can 
change the nature of the QCD phase transition [32| and 
produce gravitational waves |33[, leave its imprints on 
BBN and CMB HEl. 
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FIG. Al: The relative change of helicity (solid line) as compared to the solution of (I16|l (dashed line). The dot-dashed line 
shows the evolution of Tb in the absence of chemical potential (solely due to magnetic diffusion with lg(?7o-) ~ 16.1). A black 
dashed-dotted line shows the evolution of 7(77), growing at rf at small times and then starting to increase exponentially fast. 
A thin black vertical line marks the solution of fTS}. 



Supplementary material 



Appendix A: Evolution of a single mode 

This Appendix provides numerical results, illustrating analytic results (fT4|) . (fT6|) . and (fT8|) . Fig. I All shows the 
relative change of Tb from its initial value and the evolution of 7(77). We see that for 7 <C 1 the value of Tb (i-e. 
the total energy density of the CMF) decays over a much longer time, r]^/"/, than that of magnetic diffusion, and, a 
significantly positive value of the chemical potential (which would decay in the absence of magnetic field in a time 
77 ~ rj 1 ) is maintained during this time. 

The evolution of the mode is therefore similar to the one discussed above, as long as 7 1. For 7 ~ 1 the evolution 
quickly becomes that of standard MHD without a difference of chemical potentials, with Tb decaying due to magnetic 
diffusion and Afi dissipating through perturbativc flipping. 



Appendix B: Evolution of several discrete modes (Tf — 0) 

1. Analytic estimate of the time of the energy /helicity transfer between two modes 

The results for the time of transfer of the energy and helicity between two modes can be derived in a different way 
than done in the main text. Namely, consider the initial spectrum 

H(r] ) = H°S(k - hi) + H%5(k - k 2 ) (Bl) 

The spectrum preserves its shape throughout evolution and the tracking solution of Eq. (|19p is given by 

, . 2tt k?Hi + k%H 2 , s 

A/itr *7 « — h J r\ h T r B2) 

In the presence of such a A/x, the two modes in (jBlj) evolve as follows: 

2/«i k 2 k-\ — k 2 ■ ■ ,„ , 

H 1 = — , „ , , \r HiH 2 , H 2 = -H X B3 

a c kiHi + k 2 H 2 

To see how fast the helicity (energy) gets transferred from the mode k\ to k 2 let us combine two equations (|B3j) into 
a single equation for the ratio h(rj) = Hi (rj) / H 2 (77) (such that h(rjo) = hi > 1 and h(ij) — > as t — > 00): 

*._2a^_„ifi±« ( B4) 

at a h + q 
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(where q = fe/fei < 1)- As a result we can obtain the time difference for h to change from hi to h: 

a 



At = 



2k\ 



1 hi , 1. 1 + 
log — + - log — — 

1 — q h q 1 + h 



We can see that the time when helicities of both modes equalize (i.e. h(rj) = 1) is given by 

a 



At H = 



log hi 1 1 + /i 4 

+ - log — — 

1 -q q 2 



(B5) 



(B6) 



From Eq. (|B5P one can also easily determine the time when the energies of two modes equalize (i.e. when h = i). 

In the case when two modes with very different wave numbers (i.e. when q <C 1), the time Atu is approximately 
equal to 



a 1 

2kJ q 



2na 
kik- 



log h i 



(B7) 



Notice, that this time is much longer than the magnetic diffusion time ^fr and it depends on the ratio of initial 
amplitudes only logarithmically. 



2. Approximations, used in derivation of Eq. (1201) 

In this Section we provide some details regarding the derivation, used in the main paper and leading to Eq. (|20j) . 
Consider the system of Eqs. (flT))) . Let us assume that kiHi(r]o) 7S> ^2-^2(^0) and k 2 < k%. Then, we can rewrite the 
tracking solution (|B2[) as 

A Mtr « tatL (i _ e ) (B8) 
a 

(where e = < 1). Its time evolution is given by 

dA » = -r B (A,- 2 -^(i - <V ))) (B9) 



dt \ a 

The evolution of two modes Hi and H% is then given by 



and 



Hl = -fMjj, = -J^H 2 (BIO) 



H 2 ( ki(l-e(ri)) 1 

from which it follows that 



ff2 = d ^ir^- 1 <B11) 



fa(T?) = H 2 ( Vo )e^— (B12) 

This result is confirmed by numerical solution in Fig. lB2| where the solution of Eq. ()B12p is shown in black dot-dashed 
line (see also IBlj) . From here we can find the equation for Tb- 

2c A a*\ fa^fa fafeHa (B13) 



3. Evolution of several discrete modes 



The evolution of 10 modes, sampling a continuous spectrum with oc k 3 without perturbative chirality flip 
(solid lines in Fig. IB1[) and its comparison with the evolution of the same spectrum but for two modes only with the 
same helicities initially stored in each mode - short-dashed lines in Fig. IB11 Fig. IB2I compares the exact numerical 
solution of the case with two modes only (red and blue lines) with the result ([2"U|) of the approximate (in TLk 2 /'Hk 1 ) 
approximation (black dot-dashed line). 
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FIG. Bl: Comparison between evolution of chemical potential (left) and helicity modes (right) for 10 modes (solid lines) and 
2 modes (1st and 5th of the previous set) - dashed lines. 
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FIG. B2: Comparison of the evolution of a long-wavelength mode (in the 2 mode solution shown in Fig |Bl[l with its linear 
approximation, given by Eq. (|20[) - dashed-dotted black line. 



Appendix C: Inverse cascade and the shape of the resulting spectrum 



In the case of the continuous spectrum while 7 1 the evolution is identical to the case of Tf = 0, described above 
(see Fig. [5]). The condition 7(77) < 1 determines the time over which the described effect is present in plasma. Wc 
see from Eq. (|T7)) that for re > the fields, generated very early (probably at To ~ 100 GeV) can survive until 

T ~ 100 MeV and during the same time the non-zero chemical potential remains in plasma. The value of k at which 
the evolution stops can be determined via (|18[) with 7(7?) ~ 1. We find that k m k max (rj) - the wave- number at which 
the spectrum is peaked at the moment r\ is of the order of final time, determined by ([T7| . i.e. J7 .(fc max (r;)) ~ r\. Thus 
the value of fc max is approximately the same as in the standard MHD case, when all modes with r) a (k) < t would be 
erased by the magnetic diffusion. However, the total energy of the spectrum (or, equivalently Hfc max is much higher 
in our case - the presence of chemical potential allows for an "inverse cascade" process - the energy stored in short 
wave-lengths modes to be transferred to the longer wave-lengths (as Fig. IC1|) . 
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FIG. CI: Left: Initial helicity spectrum T-L% oc k (black dotted line), the evolved spectrum Hk at T ~ 150 MeV when 7 = 1 (red 
solid line) as well as the helicity spectrum evolved solely due to the magnetic diffusion (green dashed line) for tb ~ 5 x 10 _J . 
Right: The relative change of the total energy with (red, solid) and without (green, dashed) the chemical potential difference. 
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FIG. Dl: Ratio of chirality-flipping rates to the Hubble rate for T in GeV range (blue solid line). Red dashed line is the flipping 
rate due to weak reactions, and green (dotted) is the rate due to electromagnetic processes. 



Appendix D: Perturbative chirality flipping rates 



Fig. ID1I shows the ratio of perturbative chirality-flipping rates due to weak or electromagnetic reactions to the 
Hubble expansion rate (as a function of time) . 



Appendix E: Derivation of equations for Hk 

In this Appendix we provide the details of derivation of Eqs. <j6j) — (jTj) . 

Due to 3d translation and rotation invariance, we define the Fourier modes B k and A k for the magnetic field and 
its gauge potential in the usual way 



B(x) 



d 3 k 



ik-x r> 

e &k 



and introduce magnetic helicity density Hk and magnetic energy density p k in the fc-spacc: 

f d 3 k -> [ 

u[j]) = J j2^ Ak ' B * k = J dkUk 



(El) 



(E2) 
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FIG. Fl: Left: Evolution of Ap for two spectra with the same total energy density rs ~ 5 x 10 -5 , the same initial spectral 
index n a — 3, and different range of initial k. lg^o- = 18.7 (green), lg^o- = 16.7 (purple). Right: Zoom on initial region 
of evolution of Ap for two spectra with the same total energy density tb ~ 1 and the same initial spectral index (n s = 5) but 
different range of initial k over which this energy is distributed. 



1 f d 3 k -2 f 
PB ^2Vj (2^l Sfc l = J dk Pk (E3) 

( Bl = B-k as B is real). Notice that definitions of Hk and pk in Eqs. (|E2j) -(|E3 j) contain integrals over the absolute 
value of the 3- vector k only (cf. fl8l - l20l |. Multiplying Fourier version of Eq. ([5| by the complex-conjugated mode B%, 
we get: 

S k Bl + hlB k = - — BlB k + i — ^Bl ■ (k x B k ) (E4) 



With the use of Bf. = ik x Ak we obtain from Eq. (|E4 



If fields are maximally helical, i.e. pk = |"Hfc, Eq. (|E5j) reduces to 



dpk 2fc 2 a A/i 2 

^— = Pk + k U k (E5) 

or) o c Zir o c 

dU k 2k 2 2a Ap 

= H fc H p fc (E6) 



Appendix F: Dependence on the parameters of the initial spectrum 

Figs. [F]L show the dependence of the chemical potential difference on the range of modes k in the spectrum. 



Appendix G: Notations 



1. Magnetic energy density 



In this Appendix we discuss several conventions of expressing the energy density of the magnetic field ps = 
2^ / d 3 xB 2 , used in the literature. 
One possibility (used in this paper) is to express it in terms of the total radiation energy density 



g S .T« (Gl) 
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Alternatively, one can express ps hi terms of the total entropy of radiation s 

= P]L 

T 9 ~ ^4/3 

To convert between the tb and r g one can use: 

~ "> -no AOON"" 

2. Dimensionless quantities 

Eqs. ([4}l5l) do not change in the expanding universe with FRW metric: ds 2 = —dt 2 + a 2 (t)dx 2 — a 2 (ij)(~dr] 2 + dx 2 ), 
provided that we work with conformal coordinates x, use (dimensionless) conformal time r\ = (where we have used 
a(t) = 1/T and neglected that c/* changes with time) and substitute Afi,a, T, Tf with their conformal counterparts 
(A/ia), (era), (Ta), (F/a V a nd introduce conformal electro-magnetic fields E — > a 2 E = E c , B — > a 2 B = B c . Notice 
that aa = a c ~ const [2l|. We use these coordinates throughout the paper, starting from Eq. ((5J. To restore the 
dimensions it is sufficient to change any quantity as follows: Ap — > Afi/T, k —> k/T, Hk — > Hk/T 3 , pk —> Pk/T 4 , etc. 



= 2 W9* t3 as follows (cf. pH): 

(G2) 
(G3) 



